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We present a protocol to simulate the correlations implied by nonmaximally 
entangled two qubit states. We extend this protocol to simulate the non-local part 
of these correlations. These protocols use single cbit communication and a single 
use of Millionaire box (M-box). To the best of our knowledge, these resources are 
weaker than those used in previous protocols using classical communication. 
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I. INTRODUCTION 

One of the most intriguing features of quantum physics is the non-locality of correla- 
■^ . tions obtained by measuring entangled particles. These correlations are nonlocal because 



they are neither caused by an exchange of a signal, as any hypothetical signal should travel 
faster than light, nor are they due to any pre-determined agreement (shared randomness) 



O I as they break Bell's inequalities [1\. 

5!? ' A natural way to understand these correlations is to classically simulate them ( known 

as simulation of entanglement) using minimal resources. Obviously, this cannot be done 
using only local resources, that is, using shared randomness. The local resources must be 

rS I supplemented by non-local ones. A simple non-local resource is communication of infor- 

c^ I mation via classical bits (cbits), which we can quantify and thus provides a 'measure of 

non-locality'. 

In this scenario, Alice and Bob try and output a and (3 respectively, through a classical 
protocol, with the same probability distribution as if they shared the bipartite entangled 
system and each measured his or her part of the system according to a given random 
Von Neumann measurement. As we have mentioned above, such a protocol must involve 
communication between Alice and Bob, who generally share finite or infinite number of 
random variables. The amount of communication is quantified |[2] either as the average 
number of cbits C{P) over the directions along which the spin components are mea- 
sured (average or expected communication) or the worst case communication, which is 
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the maximum amount of communication C^(P) exchanged between Ahce and Bob in 
any particular execution of the protocol. The third method is asymptotic communication 
i.e., the limit limn^^C{P^) where P"" is the probability distribution obtained when n 
runs of the protocol carried out in parallel i.e., when the parties receive n inputs and 
produce n outputs in one go. In this paper we are concerned with the worst case sce- 
nario. A fundamental result for this scenario is that kT^ ( k a constant) cbits of classical 
communication is required to simulate the correlations implied by a ra qubit maximally 
entangled state [3]. This was followed by a remarkable result due to Toner and Bacon [1] 
who showed that a single cbit of communication is enough ( apart from shared random 
variables) to simulate the correlations of a two qubit singlet state. We have later shown 
[5], [6] that for simulating entanglement of arbitrary spin S singlet state, communication 
of n = log2(s + 1) cbits is enough, provided we confine only to spin measurements ( of 
type a ■ S). 



II. SIMULATION USING NON-LOCAL BOXES 

Another fruitful approach to this problem is to use PR-box [^. PR-box is a concep- 
tual and mathematical tool developed to study non-locality, first proposed by Popescu 
and Rohrlich [8j. It was demonstrated in [9J, that the correlations of the two qubit singlet 
can be simulated by supplementing hidden variables ( shared randomness) with a single 
use of the PR-box. Although mathematically based on the Toner and Bacon [Ij result, 
this work is a major conceptual improvement, as the PR-box is a strictly weaker resource 
than a bit of communication, because it does not allow signaling. 

However, it was soon discovered that a single use of PR-box is provably not sufficient to 
simulate some partially entangled two qubit states [10]. Simulation of these states has 
been one of the open problems of the whole field. For example, [11] , [12] use a signaling 
resource strictly stronger than one cbit of communication. A reason for this apparently 
surprising result is that PR-boxes have random marginals, ((a) = = (/5), a fact which 
is consistent with two qubit singlet state) while the correlations arising from partially 
entangled quantum states have nontrivial marginals. 

Thus it appears that it is especially difficult to simulate at the same time nonlocal corre- 
lations and nontrivial marginals, like these corresponding to partially entangled quantum 
states. Recently, Brunner, Gisin, Popescu and Scarani [13] (hereafter referred to as BGPS) 
have given a procedure to simulate entanglement in non-maximally entangled states us- 
ing four PR-boxes and one M-box (see bellow). In order to overcome the above difficulty, 
they introduce the concept of correlated local flips, which is independent of whether we 
use nonlocal boxes or classical communication to simulate entanglement. The idea is 
that some nonlocal box, or a classical communication protocol, flrst simulates non-local 
correlations with trivial marginals and then use local flips to bias the marginals. In the 
following we outline the simulation scheme, based on local flips. 
The problem is to simulate the quantum correlation implied by a general non-maximally 



entangled two-qubit state 

|z^(7)) = cos(7)|00) + sin(7)|ll); (0 < 7 < \). (1) 

Alice and Bob perform measurements along directions a and h on their qubits ( that is, 
they measure operators S ■ a and S ■ b). Let a, P {a, P E { — 1, +1}) denote their outputs 
respectively. The joint probability Pqm{o:, P\d,b) is given by 

Pqm{o:, P\a, b) = -[1 + acttz + Pcb^ + ajSC^a, 6)] , (2) 

where c = cos 27 and s = sin 27 and 

C{d, b) = {ap)QM = {azbz + s{a^b^ - ayby)) (3) 

In order to simulate this joint probability, we proceed as follows. We first set up a 
procedure ( involving non-local boxes or communication) to simulate the joint probability 
PQ{a^(3\a,b) = | [l + a/3Co] where Cq is the correlation {ajS) = '^aPPoi^a, P\a,b) before 
flipping. 

Now, following BGPS, we invoke the local flip operation as follows. Alice (Bob) flips the 
output -1 with probability fa{fb) while the output +1 is left untouched. After the local 
flipping operation with probabilities fa and ft, respectively by Alice and Bob, assuming 
fb^ fa, the joint probability Po becomes 

Pf{a, P\d, b) = ^[l + afa + Pfb + aPUa + (1 - fb)C,)\ . (4) 

In order that Pf coincides with Pqm-, we identify /„ = cuz and fb = cbz- Note that 
the condition fl > fa now becomes bz > az- This gives, Cq = a ■ B where B = 
{sb^, sby, bz-c) / (l- cbz). If fa > ft, {a-z > bz) , then 

Pj{a,P\a,b) = ^[1 + afa + Pfb + aP{ff,+ {I - fa)Co)]. (5) 

where Cq = A ■ b and A = {sax, say, a^ — c) / (1 — caz). It is easy to see that A and B 
are unit vectors. The joint probability Pq to be simulated before local flipping operation 
depends on whether {bz > az) or (a^ > bz) { via Co). But Alice (Bob) cannot have any in- 
formation on bz {az). In order to pave way through this situation BGPS invoke a nonlocal 
box called M-box which has two real inputs x,y E [0, 1] and binary outputs m,n E {0, 1}. 
The M-box is defined hj m(Bn = [x < y] where [x < y] is the truth value of the predicate 
X < y for given values of x and y. 

III. THE PROTOCOL 

In this paper we give a protocol for simulating entanglement in an arbitrary non- 
maximally entangled two qubit quantum state, which uses single cbit of communication 



and a single use of M-box. The protocol runs as follows. 

Alice and Bob share seven random unit vectors, distributed independently and uniformly 

over unit sphere, Ai, A2, /ii , /i2, As, /i4, /^s- Also Alice and Bob share a common reference 

direction z. 

(i) Alice and Bob input a^ and h^ values respectively in M-box, which outputs m and n 

as above. We replace m and n by p and q respectively as p = 2m — 1 and g = 2n — 1 so 

that p, gG{— 1,1}. Alice gets p without knowing q and Bob gets q without knowing p. 

(ii) Alice outputs 



a. 



sgn[ux-\\\^\ —r— 1 sgn ( M2 ■ Ai 



(6) 



where 111^2 = sgn {z ■ jli^^) a + sgn {z ■ jl2,z) -^ + a-i,2- Here 01^2 are vectors to make Ui^2 unit 

vectors. 

(iii) Alice sends one cbit to Bob 

Ca = I -^:r- I sgn (ui ■ Ai ) sgn (ui ■ A2 ) + ( -—— ) sgn ( M2 ■ Ai ) sgn ( U2 ■ A2 ) (7) 



(iv) Bob outputs 






vi- Ai + CaA^ 



+ I —^- I sgn 



^^2 ■ Ai + CaA2 



where 'Oi_2 = sgn {z ■ (13^2) b+sgn {z ■ fii^^) JB+sgn {z ■ ft^) bi^2- Here 61^2 are vectors to make 
■1)1^2 unit vectors. Since / sgn (mi,2 ■ Ai ) dXidX2 = = / sgn 'Oi_2 ■ ( Ai + CaA2 ) dXidX2, 
we see that at this stage of the protocol (a)Q = = (/3)o- For {aP)^ we have, after 
integrating over A 1^2 



(«/?) 



{A} 



1+P + q + pq 


4 
1+p-q-pq 



Ui- Vl + 



1 — p — q + pq 

4 



U2 ■ f 2 



, , . . . ... .. , , l-p + q-pq ... .„. 

+ I : I Ml ■ ^2 + ( ]U2-Vi (9) 



Next, integrating over {ft} we get 






(10) 



(v) Alice and Bob perform the local flip operation with probabilities fa = cttz and fb = cbz 
respectively. Now, if /;, > fa{bz > dz) then p = q which implies from Eg. (1101) that 
(a/3)o = a ■ B = Cq, so that (a/3)^ = ca^ + (1 - cbz)d ■ B = {af3)Qjy.j. If /« > /b(a^ > b^) 

then p 7^ g implying (a/3)o = A ■ b = Cq, so that (a/?) j = cb^ + (1 — caz)A • b = (a/?)^^^/- 



The local flipping operation ensures that {a) r = fa = caz = {c()qm ^^^ (/?)r = fb = 

chz = {13) QM- We see that the protocol simulates the joint probability Pqu^o., [3\a,h) as 
in Eq.©. 

IV. EPR2 DECOMPOSITION AND SIMULATION OF NON-LOCAL COR- 
RELATION 

Elitzur, Popescu and Rohrlich (EPR2)[Tlj considered an experiment involving many 
photon pairs. The question they address is whether a subset of these pairs can be de- 
scribed with local correlations only, while the remaining ones are described non-locally. 
The global statistics, comprising the local and the non-local subsets of pairs, should repro- 
duce the quantum statistics. Formally the EPR2 approach is to decompose the quantum 
correlations Pq as a convex sum of a local probability distribution P^ and a non-local one 
Pnl '■ 

Pq = pl{p)Pl + (1 - Pl{p))Pnl. (11) 

The weight Pl{p) is a measure of locality of the state p. A decomposition of the form 
(ITTl) is particularly well-suited for the task of simulating quantum correlations, since 
only the non-local part Pnl has to be simulated, the local part requiring only shared 
randomness. We use the improved version of the original EPR2 decomposition presented 
in [15]. For the states \ip{'^)) in Eq.([l]) we are interested in, the weight of the non-local 
part PNiil) = (1 ~ Pl(7)) vanishes in the limit 7 — ;► of separable states. We present 
a simulation protocol for this Pnl with one use of M-box and one cbit communication 
followed by correlated local flips. This second model, involving Pl,Pnl decomposition, 
fulfills the desideratum that very weakly entangled states can be simulated by a vanishing 
amount of non-local resources. The exact form of the decomposition described above, as 
given in [i5\ is 

^W + aF{az)+pF{bz) + a(3G{a,b)\ (12) 



where 
with 



F{x) = - ex — {1 — s)f{x) 



and 



f{x) = sgn{x)min I 1, |x| 

G{d, b) = aj)^ - ayby + - [a^b^ - (1 - s)f{az)f{bz)] (13) 

Scarani [15J has shown that if both Alice's and Bob's measurement settings lie in a slice 
around the equator of the Bloch sphere, defined by 

|a.|,|&.| <^^, (14) 



Pnl takes the simple form -F(a^) = = F(h 



Pnl-^ 



1 + af3{a^bx - ayhy - a^h-^) 



1 + al3{a ■ h') 



with h' = {bj;, —by, —bz). The local marginals are random and the correlations reduce to a 
simple scalar product, as in the singlet case. Further, if only one of the conditions flT^ 
is true, then the corresponding marginal is random and the local flip operation becomes 
redundant. 

Note that, since Pqui—Oi, [3\ — a,b) = PQM{<y, P\a,b) it is sufficient to consider the case 
where az,bz > 0. Further, for each of the four possibilities corresponding to inequality 
(fT4l) it is straightforward to check that {aP)p = G{d, b) where {aP)p is the correlation 
after flipping operation. We now give a protocol using one use of M-box and one cbit of 
communication to simulate P/vl of Eq. (TT2!) in all the four cases, namely. 



a-z < ,bz < 

c c 

l-s 



a, < , bz > 



a^ > 



6. < 



c 



Ctz > 



c 
l-s 



c 



,bz > 



c 

l-s 
c 



(15) 



Alice and Bob share 9 random unit vectors, distributed independently and uniformly over 

unit sphere, Ai, A2, fti , /i2, As, ^4, As, Ae, At- 

i) As in the first protocol Alice and Bob input az and bz values respectively in M-box, 

whose outputs p and q are received by Alice and Bob respectively. 

(ii) Alice outputs 
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sgn ( Ml ■ Ai ) + ( -—— 1 sgn ( M2 ■ Ai 



sgn Mo- A] 



(16) 



where 



Mo = [sgn {z ■ fii) + sgn {z ■ (I4) + sgn {z ■ ftg)] a + ao. 



Here ao is vector to make uq unit vector and M1.2 are the same as those in the first protocol 
except that A is redefined as A = {sQx, say, c — az)/{l — caz). 
(iii) Alice sends one cbit to Bob 
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sgn ( Ml ■ Ai ) sgn ( Wi ■ A2 ) + ( — - — 1 sgn ( ^2 ■ Ai ) sgn [u2- ^2 



sgn [uo- Xi) sgniuo- X 



(17) 



iv) Bob outputs 
1-s 



/3 = 
+ 
where 



6. < 



c 
l-s 



1 + q 



sgn 



^^1 • Ai + c^A^ 



+ I -^^ I sgn 



^^2 ■ Ai + c^A^ 



sgn 



^^0 ■ Ai + CaX2 



vo = [sgn [z ■ fi2) + sgn (z ■ ft^) + sgn {z ■ /ie)] b' + sgn [z ■ fij) bo. 



Here bo is vector to make -Oq unit vector and -01^2 are as in the first protocol except that b 

is replaced by b' defined above. 

After integrating over all A and fi we get 



(a/5)c 





r i-si 




L 1-^1 
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r i-si 
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az < 




bz < 


a-b' 
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-a.B+^—MA.l 



az > 



l-s 



&. < 



l-s 



A-b' 



(19) 



(v) Alice (Bob) performs local flip operation with probability F{az) {F{bz))- Note that 
F(az) = ( F(bz) = 0) whenever a^ < ^— ^ (6^ < ^— ^) , so in this case, Ahce ( Bob) does 
not flip. It is straightforward to show from Eq. (TT^ that G{d, b) (Eq. (IT5]) ) equals {a(3)p in 
all the four cases [T51 Since the local flipping with probabilities F{az) = and F{bz) = 
produce required marginals, we see that our protocol simulates P^l ( Eq. flT2l) ). 

V. CONCLUSION 

In conclusion, we have presented protocols to simulate partially entangled two qubit 
states, using one cbit of communication and one use of M-box. Our work using EPR2 
decomposition Eq. OllI) naturally conflrms that the less the quantum state is entangled, 
less frequently one needs to use non-local resources as the simulation by BGPS shows. 
Further, the amount of non-local resources needed to simulate a partially entangled state 
is larger than that required to simulate maximally entangled states. The simulation of 
partial entanglement in BGPS uses all non-signaling resources [ PR and M-boxes] while 
our protocols use both the signaling and non-signaling resources. At any rate, our protocol 
is an improvement on previous protocols which used more than one cbit communication 
|llj . as our protocol uses one cbit communication and one M-box which is a weaker re- 
source. Since the existence of no signaling protocol and more than one cbit protocol does 
not imply the existence of a protocol with one cbit communication and one M-box, flnding 
one such protocol is a new discovery in its own right and hopefully will flnd applications 
in future. 
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